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Abstract 

It is well known that short-time expansions of heat kernels correlate to for- 
mal high-frequency expansions of spectral densities. It is also well known 
that the latter expansions are generally not literally true beyond the first 
term. However, the terms in the heat-kernel expansion correspond rigorously 
to quantities called Riesz means of the spectral expansion, which damp out 
oscillations in the spectral density at high frequencies by dint of performing 
an average over the density at all lower frequencies. In general, a change of 
variables leads to new Riesz means that contain different information from the 
old ones. In particular, for the standard second-order elliptic operators, Riesz 
means with respect to the square root of the spectral parameter correspond to 
terms in the asymptotics of elliptic and hyperbolic Green functions associated 
with the operator, and these quantities contain "nonlocal" information not 
contained in the usual Riesz means and their correlates in the heat kernel. 
Here the relationship between these two sets of Riesz means is worked out 
in detail; this involves just classical one-dimensional analysis and calculation, 
with no substantive input from spectral theory or quantum field theory. This 
work provides a general framework for calculations that are often carried out 
piecemeal (and without precise understanding of their rigorous meaning) in 
the physics literature. 



E-mail: fulling@math.tamu.edu 
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I. PHYSICAL MOTIVATION 



Let if be a positive, self-adjoint, elliptic, second-order partial differential operator. For 
temporary expository simplicity, assume that H has purely discrete spectrum, with eigen- 
values A n and normalized eigenfunctions <f) n (x). Quantum field theorists, especially those 
working in curved space, are accustomed to calculating (1) the heat kernel, 

oo 

K(t,x,y) = J2Mx)MyYe- tK ; (1) 

n=l 

(2) the Wightman function, 

OO 1 

W(t, x,y) = Y J Mx)MvYirn= e ~* y/ *> ^ 

n=l Z V *n 

which determines the vacuum energy density of a quantized field in a static space-time 
background. Of greatest interest are the behaviors of these functions with y = x (either 
evaluated pointwise, or integrated over the whole spatial manifold) in the limit t — *■ 0. 

Roughly speaking, the small-t asymptotics of both ([!]) and (|2]) are determined by the 
large-A (ultraviolet) asymptotics of the density of eigenvalues and of the eigenfunctions. 
However, there is a major difference between the heat kernel and the Wightman function in 
this regard. All the coefficients b s (x) in the heat-kernel expansion 

oo 

K(t,x,x)~Y.Ux)t^ H (3) 

s=0 

(m = dimension) are locally determined by the coefficient functions in the differential oper- 
ator H at the point x. The same is true of the leading, singular terms in an expansion of 
W(t,x,y), which are removed by renormalization; the finite residue, however, is a nonlocal 
functional of the coefficients in H, containing information about, for example, the global 
structure of the manifold — this is what makes vacuum energy interesting and nontrivial. It 
follows that this nonlocal information is somehow contained in the ultraviolet asymptotics 
of the spectrum, although it is lost in passing (pointwise) to K(t,x,x). The Wightman 
function is more typical of Green functions associated with H; the striking and somewhat 
mysterious thing is how special the heat kernel is. 

The primary purpose of this paper is to point out that many of the facts of this subject 
have nothing to do specifically with partial differential operators (much less with quantum 
field theory). Rather, they result from some classical real analysis (in one dimension) con- 
cerning the summability of series and integrals, much of which was developed by M. Riesz 
and G. H. Hardy early in the twentieth century. One consequence is that many formal rela- 
tionships between the asymptotic expansions of Green functions and those of the associated 
spectral measures, and between the asymptotic expansion of one Green function and that of 
another, can be worked out in the abstract without reference to the detailed spectral theory 
of any particular operator. The same is true, qualitatively, of the limitations of such rela- 
tionships: The rigorous asymptotic expansion of a heat kernel suggests a formal asymptotic 
expansion of an associated spectral measure, but the latter is usually valid only in some 
averaged sense. It can be translated into rigorous statements about the Riesz means of the 
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measure. The construction of Riesz means not only washes out oscillations in the measure 
at the ultraviolet end, but also incorporates some information about the measure at small or 
moderate values of A. The information contained in a Riesz mean depends on the variable 
of integration; for example, the Riesz means with respect to v^A contain more information 
than those with respect to A. The difference between the nonlocal asymptotics of the Wight- 
man function and the local asymptotics of the heat kernel is, at root, an example of this 
phenomenon. 

These relationships can be made yet more precise by introducing some concepts from 
distribution theoryoH 



II. SURVEY OF THE PHENOMENA 

To provide a concrete context for the later sections of the paper, we display here some 
Green functions and asymptotic expansions associated with a constant-coefficient differential 
operator in one dimension. (The later sections are logically independent of this one and much 
more general). 



Here and in Sections |V| and [VII| we use as a surrogate for the Wightman function a 
technically simpler Green function that we call the cylinder kernel. This is the integral 
kernel T(t, x, y) of the operator e _i ^, and it is related to the elliptic equation 

■w = ™ < 4) 

in the same way that the heat kernel (the integral kernel of e~ tH ) is related to the parabolic 
equation 

-- = **. (5) 



That is, 



*{t,x)= I T(t,x,y)f(y)dy (6) 

J M 

solves (|J) in the infinite half-cylinder (0, oo) x M. with the initial data 

lim*(t l x) = /(x) (7) 

on the manifold M.. The cylinder kernel displays the same kind of nonlocal short-time 
asymptotics as the Wightman function. 
In this section we consider 

d 2 

on various one-dimensional manifolds M.. 
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Case M = R. 



The heat kernel is 



As 1 1 0, 



K(t,x,y) 



K(t,x,y) = (A7it)- 1/2 e~ ix - y)2/u . 

\ (4vrt)- 1 /2 + (t°°) ify = x, 



0{t c 



ify^x; 



that is, all terms in the expansion (|3[) beyond the first vanish. 
The cylinder kernel is 



t 1 

T(t,x,y) = — 



7r (x — y) 2 + t 2 



As t[ 0, 



T(t,x,y) 



1 

7Tf 



if y = z, 



E(-i) fc/ 



2A- 



k tt(z - y) 2 ^ ' 

(For the distributional, rather than pointwise, limit, see Ref. [I].) 



if y 7^ x. 



(9) 



(10) 



(11) 



(12) 



Case M = H + 

We consider the operator (H) on the interval (0, oo), with either the Dirichlet or the 
Neumann boundary condition at 0. The Green functions are most easily obtained by the 
method of images from the previous case. The heat kernel is 



K(t,x,y) = (47it)- 1 / 2 



e -(x-y) 2 /4t e -(x+yf/At 



(13) 



where the upper and lower signs are for the Dirichlet and Neumann cases, respectively. 
Because of the rapid decay of the image term, the asymptotic behavior is still described 
exactly by fllPf) ; the heat kernel in the interior does not sense the existence of the boundary. 
The cylinder kernel is 



T(t,x,y) 



7T 



=F 



1 



x-y) 2 + t 2 (x + y) 2 + t 2 



(14) 



As t[ 0, 



T(t,x,y) 



nt T n(2x) 2 ^ [ ' \2x) 



2k 



iiy = x, 



2k 



ir(x-y) 1 



/c=0 



X 



=F 



yj n(x + y)^ 



t 00 ft 

' E(-i) fc 



2/v 



z + y 



if y ^ x. 



(15) 
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Because of the slow decay of the basic kernel in its role as image term, the expansion differs 
from fll2|) beyond the leading 0(t _1 ) term; the cylinder kernel senses the presence of the 
boundary, the type of boundary condition, and the distance to the boundary (more precisely, 
the length, x + y, of a path between the two arguments that bounces off the endpointi). 



Case M = S 1 

Consider (|8|) on the interval LrL, L) with periodic boundary conditions. The heat kernel 
is the well known theta functionQ 



K(t,x,y) = (4nt)- 1 / 2 ]T e 



-(x-y-2nL) 2 /At 



(16) 



The expansion (PH| ) is still valid. The cylinder kernel could also be expressed as an infinite 
image sum, but its Fourier representation can be expressed in closed form (via the geometric 
series) : 



T(t,x,y) = ^1 ( 

k=—oo 



1 00 

1 V J-*k(x-y)/L -K\k\t/L 



sinh(7r£/L) 



2L cosh(7rt/L) — cos(7r(x — y)/L) 
The first few terms of the expansion as t [ are 



(17) 



T(t,x,y) 



r 1 

Tit 1 



1 + 12(l) -™>(t) +0(() 



[ 2L 2 1 — k 



1 + 



Tit 



LJ \6 2(1 — k) 



+ 0(t 4 



if y = x, 
iiy^x, 



(18) 



where 



k = cos • 



n( x - y) 

L 



(19) 



Thus the cylinder kernel is locally sensitive to the size of the manifold. (In the limit of large 
L, (|l^) matches (0), as it should.) 

In this case it is possible to "trace" the Green functions over the manifold: 



K(t) = J L ^ K(t, x, x) dx ~ (4vrt)- 1/2 (2L) + 0(t 

(since ( |TUD is uniform in A4), and 

m / x f L , 1 sinh(7ri/L) 2L ' 

T{t) = J L T{t,x,x)dx 



(20) 



cosh(7rt/L) — 1 TTt 



1 /7lt\ 2 



(21) 



Both have leading terms proportional to the volume of the manifold, but (|2l|) has higher- 
order correction terms analogous to those in fljjl) . 
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Case M = I 



We consider (|5p on the interval / = (0,L). For brevity we consider only Dirichlet 
boundary conditions and the expansions on the diagonal (coincidence limit). The heat 
kernel is 



CXJ 



-y-2nL) 2 jit _ e -(x+y-2nL) 2 /it 



(22) 



Expansion fllQf) is valid in the interior (but not uniformly near the endpoints). The cylinder 
kernel is 



T(t,x,y) 



2L 



sinh(7rt/L) 



sinh(7ri/L) 



cosh(7rt/L) — cos(7r(x — y)/L) cosh(7rf/L) — cos(7r(x + y)/L) 



Not surprisingly, its expansion combines the features of fll5|) and (|i~8|). 



T(t, x, x) 



7rt 



1 /?rt 
1 + l2lT 

7rf 1 



1 /7rt 



2L 2 1 - cos(2ttx/L) 



720 V L J 
1 



+ 



id 



Lj I 6 2(1 -cos(27nr/L)) 



+ ■ ■ 



1 

+ 



LJ V 12 2(1 - cos(27ra;/L)) 



7Tt\V-l 1 



+ 



L7 V720 12(1 - cos(2ttx/L)) 4(1 - cos(2vra;/L)) 2 / 
(Compare ( |i5|) in the form 



+ 



T(t, x, x) 



id 



4x 2 + 16x 4 ^ 



which (|24|) matches as x — ► or L — > oo. 
For the traces in this case one has 



= / L i^(t, x, x) cfe ~ (4vrt)- 1/2 L - \ + 0(t c 
Jo 



a well known resultn, and 

r-L 



T(t) = / T(t, x, x) dx 
Jo 

1 sinh(7rt/L) sinh(7rf/L) r L 



dx 



2 cosh(7rt/L) — 1 2L Jo cosh(7rt/L) — cos(27ra/L) 

1 sinh(7rt/L) 1 

2 cosh(vrt/L) - 1 ~ 2 
L 

7rt ~ 2L 12 V L 



i-^-^f?y+o(t 4 ) 



(23) 



(24) 



(25) 



(26) 



(27) 
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In comparison with ( pQ|) and (pT|), the leading terms of (|26|) and fl27|) have adjusted to reflect 
the smaller size of the manifold, and the new second (t-independent) terms are the effect of 
the boundary. The minus sign on those terms distinguishes the Dirichlet boundary condition 
from the Neumann. (A less trivial differential operator would yield more complicated ex- 
pansions, with higher-order terms exhibiting an interaction between the boundary condition 
and the coefficients (potential, curvature) (see Ref. ||).) 

It is worth noting that the most important dependence of ([T!]) or ( [2"7j ) on the size of 



the manifold comes from the integration, not from the form of the integrand. (Indeed, 
the L-dependence of flT8"| ) and fl2"l| ) as written is downright misleading in this respect.) In 
ignorance of both the size of Ai and the nature and location of any boundaries, knowledge 
of the 0(t) term in T(t, x, x) on a small interval of x would be a rather useless tool for the 
inverse problem. 

III. NOTATION AND BASIC FORMULAS 

The basic references for this section and much of the next are HardyiHl and Horman- 
dei@0. (The formulations given here are somewhat new, however.) 

Let /(A) be a function of locally bounded variation on [0, oo) such that 

/(0) = 0. (28) 

Typically, / will be defined Stieltjes integral 

/(A) = [ X a(a)dn{a), (29) 



o 

where /i(A) is another function of the same kind, and a(A) is (say) a continuous function. 
By convention we take functions of locally bounded variation to be continuous from the left: 

/(A)=lim/(A-£); (30) 

rb rb—e 

I a(a) dfx(a) = lim / a(<r) d/i(a"). (31) 

By &\f we denote the cvth derivative of / with respect to A. Since iterated indefinite 
integrals can be regarded as derivatives of negative order, we define 

d^ a f(X) = fdaf- do a j{a a ). (32) 

Jo Jo 

As is well known, the iterated integral is equal to the single integral 

<9a Q /(A) = /"(A - a)- l f(a) da, (33) 

(a — 1)1 Jo 

which in turn may be converted to 

d- x a f(\) = ±j\\-ardf(a). (34) 
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This last formula remains meaningful for a = 0, yielding the natural definition 



d° x f(X)= f X df(a) = f(\). (35) 
jo 

The ath Riesz mean of / is defined by 

R<*f(\) = a\ \- a d- x a f{\) =J Q \l- j) a df(cr). (36) 

We call d^ a f(X) the ath Riesz integral of /, and we call 

lim R«f(\), (37) 

A— too 

if it exists, the ath Riesz limit of /. If such a limit exists for / defined by (|29f) , we say that 
the integral ( p9|) is summable by Riesz means of order a. Eq. ( p4[ ) may be used to define 
all these things for any a GC with Re a > —1, but here we shall consider only nonnegative 
integer a. 

From ( p2|) and (|33|) with / replaced by d^f, we have (for a > 0, /3 > 0) 



/(A) = T-^rry /"(A - a)- 1 ^/^) da, (38) 
(a — 1 ! Jo 



hence 



R«+Pf(\) = « ^ j A" (Q+/3) / o (A - ( r) Q - 1 ^^/(a) da. (39) 

From this is proved Hardy's first theorem of consistency: 

If an integral (29) is summable by Riesz means of order «o , then it is (40) 
summable to the same value of means of order a for any a > «o • In 
particular, if the integral is convergent, then it is summable by means 
of any order a > 0. 



The proof is similar to that of (|45|), below. 

Although Riesz and Hardy were primarily interested in defining numerical values for non- 
convergent series and integrals, our concern (following Hormander) is to use Riesz means as 
a way of organizing the asymptotic information contained in a function, without interference 
from small-scale fluctuations. Thus we are more interested in the Riesz means themselves 
than in the Riesz limits. 



IV. THE RELATION BETWEEN MEANS WITH RESPECT TO DIFFERENT 

VARIABLES 

Let uj be related to A by 

X = u k (k > 0, k ^ 1). (41) 
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We are primarily interested in the cases k = 2 and k = ~. (More generally, one could treat 
two variables related by any orientation-preserving diffeomorphism of [0, oo).) We write a 
and r as integration variables corresponding to A and u respectively — hence a = r k . We 
let 

/» = f(u k ) = /(A) (42) 

and may omit tildes when no confusion seems likely. 

The Riesz mean B%f can be expressed in terms of the mean R%f [ = R"f] and vice 
versa. We call the following the Hardy formula because it is implicit in Hardy's 1916 proof 
of the "second theorem of consistency" : 

Rtf = k a R a J + f" J k>a (\, a)R a J da, (43) 
J o 

where 

In (pE3"|), R®f means R^f evaluated at r = o x l k . In (Hp, the ratio of T functions is interpreted 
as if kj + k — a is a nonpositive integer. 

Proof: Use successively the definitions (|3~6"D, (|3~3"1) , (|4lD , and (j32j): 

= aA" a / (A — cr)" _1 /(cr) dcr 
Jo 

= au~ ak r{u k - r^^kT^fir) dr 
Jo 

= au~ ak f (c/ - r*) a - 1 Ar fc - 1 a?(ar a[ /) rfr. 

Jo 

Integrate by parts a times. The lower-endpoint term always vanishes, because d~^f(0) = 
for all (3. Until the last step, the upper-endpoint term contains positive powers of uo k — r k , 
so it also vanishes. At the last step there is an upper-endpoint contribution 

R x = {-ly^auj-^ia - 1)! (-kr k - 1 ) a {-l)d; a f = a\ k a u~ a dZ a f = k a R a J. 

T=LJ 

The remaining integral is 

R 2 = (-\) a auo~ ak d*[(u k - T k ) a - 1 kT k - 1 ]d~ a f dr. 
Jo 



But 



Q-l 



a — I 



oj — r kt — k > — i J , \u v J 'r 



j=0 

so 
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-l) a au- ak d*[{uj k - T k ) a - l kT k -^ 



a-l 

akJ2(-±) a ~ j; 

j=0 



a — 1 

, 3 , 



a-l 

U(kj + k-i-i) 

.i=0 



kj+k-l-a^-kj-k 



Thus, by fl36|) applied to /(r), 

fl2 = rf -±vr% T^A^-* i Kfdr 

Jo p o j\(a-l-j)\T(kj + k-a) V ; a\ rJ 

[ X J Ka {\a)R a T fda. 
Jo 



Adding i? x and R 2 , one obtains the formula (|43| ) to be proved. 
As a corollary we see: 

If an integral (29) is summable by ath-order Riesz means with respect (45) 
to A, then it is summable to the same value by ath means with respect 
to u, and conversely. 

Proof: We show the converse; the direct statement follows when one replaces A; by 1/k 
and interchanges A and u. It suffices, since Riesz means are linear, to consider the special 
cases (i) where the Riesz limit in question is 0, and (ii) where /(A) = C, a constant, for all 
A / 0. The latter case is trivial since all Riesz means and hence all Riesz limits are equal 
to C . In case (i), we are given that R°f — > as a — > oo; that is, for every e > there is a 
K such that \R°f(a)\ < e when a > K. For A > K, write 

R 2 = f K J k , a (\ 7 a)R%fda+ f X J k . a {\a)R a T fda = R 21 +R 22 . 
Jo Jk 

From the form of Jk, a (PD; R21 approaches as A — > 00, and |i?22| has the form 

\R22\ < eJ^CjX'i- 1 ( X a j da = eO(A°), 
3=0 " Jk 

which can be made arbitrarily small. Thus R%f approaches 0, as asserted. 

A generalization of this argument (due to Riesz and HardyQ'Q) proves the second theorem 
of consistency (for a integer): Let A = g(u), where g is an increasing function in C°°(R + ), 
g(0) = 0, g(oo) = 00, and 

ar u g(u) = 0(u- r g(u)) for all r = 1, 2, . . . . (46) 

Then ( ^5|) holds, except possibly for the clause "and conversely". Conditions sufficient 
together to guarantee (^) (given the other conditions on g) are 

g(cu) = 0(lu a ) for some A > (47) 

and 
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g{oj) is given (for sufficiently large uj) by an explicit, finite formula (48) 
involving the logarithmic and exponential functions, real constants, 
and (real, finite) elementary algebraic operations. 

The significance of this theorem is clearer in the following imprecise paraphrase: (1) If uj 
increases to oo, but more slowly than A, then A-summability implies tu-summability, provided 
that the increase of uj with A is sufficiently "steady" — that is, derivatives of g must not 
oscillate so as to disrupt ([46]). (2) Under those conditions, cu-summability does not ensure 
A-summability — e.g., an integral i? a -summable with respect to to may not be i? a -summable 
with respect to A ~ e w . However, if the rates of increase of A and uj differ only by a power, 
then the two types of summability are coextensive. 

If k is an integer (necessarily > 2), there is another relation between A-means and uj- 
means, which we call the Hormander formula [Ref. [U], Sec. 5]: 



R a JW 








i=i 



df(r) 



ak 

= ^2 bpR^f f° r certain numbers bp . (49) 

This is obtained from (|36|) by expanding the factor (1 — er/A) = 1 — [1 — (1 — r/uj)] k by 
the binomial theorem. In view of (fiDD, (f49"D has the converse part of (f45|) as a corollary. 
More significantly, the absence of an integral in fl4*9| ) as compared with (H3) means that the 
asymptotic (A — > oo) behavior of B%f is entirely determined by the asymptotic behavior 
of R®f, a conclusion which we shall reach independently below. We shall see also that the 
converse is false: The asymptotic behavior of R°f is affected by the values of / at small 
A, in a way that is not captured by the asymptotic behavior of R^f for any (3, no matter 



how large. Thus it is really essential in (|49"D that k is an integer (so that the binomial series 
terminates). 

We now work out in each direction the relation between the Riesz means with respect to 
the spectral parameter and those with respect to its square root, taking note of a fundamental 
difference between the two calculations. 

1. First assume that as A — > oo, 



m a 1 . 



i?>(A) = £ a as \--? + 0(A----2 ) (50) 

s=0 

for some positive integer m, as suggested by Hormander's results on the Riesz means of 
spectral functions of second-order operators (where m is the dimension of the underlying 
manifold). Let us attempt to determine the asymptotic behavior of R™jl (A = uj 2 ) from (43) 
and ( fHD with / = k — |, and A and to interchanged: 

= 2- a J2SA*+ r J ha (uJ,r)RydT, (51) 

JO 2 ' 

A>.^g(- 1 r-^-' j!( J 1 _. )!r( ^tjV - (52) 
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Note that all terms in (|5^) with j odd are equal to zero, since + 1) is an integer < a. 

Since \i is of locally bounded variation, it is bounded as A j 0. Hence d^ a ji = 0(X a ) at 
small A (see Remark at end of Sec. |V|), so is bounded there. Thus there is no problem 
with convergence at the lower limit of the integral in (|5T1), as a whole. However, some of 
the individual terms in (|50|), and hence the remainder term by itself, will be singular as 
t I if a > m. Therefore, for each integral encountered when (|50| ) is substituted into ([51~1), 
one must choose an appropriate lower limit of integration, r . Let R be the part of the 
total integral thereby omitted; thus the integrand of R equals Ji a i?"/i(r 2 ) for small r, and 
is redefined as r increases — in such a way that R is convergent and depends on u only 
through the factors uj^^ 1 . Then R will be of the form 



a-l 



R= Y, Z j^ 3 '\ (53) 



3=0 
j even 



where the Zj are constants, which cannot be determined from the information in ( |50| ) since 
they are affected by the behavior of R"n at small A. Next consider the part of the integral 



in (51) associated with the series in (50): 



Va as V (-l)"-^- 1 — r m ~ s+ ^dr. (54) 

j even 

When each integral in is evaluated, the contribution of the lower limit is of the form 
(|53|) and may be henceforth counted as a part of R. The contribution of the upper limit is 
proportional to u m ~ s , unless m — s is odd and negative, in which case the term j = — m+s — 1 
yields something proportional to uo m ~ s In u and the other values of j yield more terms to be 
absorbed into R. Finally, when integrating the remainder term in ( [50f) we may assume that 
t and us are so large that for some K, the integral is less than 



K I" \Ji ta (u,T)\r m - a - 1 dT < Y K 3 uj- j - 1 T r m - a+j - 1 dr 

j even 



a-l 



j even 

except that if m — a + j = 0, the corresponding term is 

K j w mr - a - 1 (kiu;-]nT ). (56) 

Now Tq must be chosen differently for different j: If m — a + j < 0, take r — > oo, leaving in 
(j55|) a single contribution to the error of order 0(uj m ~ a ~ 1 ). (An integral over a finite interval 
is thereby included in R, but ( |53D is still valid.) If m — a+j > 0, we have —j — 1 < m — a — 1, 
so that both terms in (|5^) are 0(u; m_a_1 ) for finite To , with a possible extra logarithmic 
factor in the worst case flSHp. 

Therefore, adding all contributions in fl5lD, one has that whenever R®H has the asymp- 
totic behavior (0), has the behavior (as uj — ► oo) 
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where 



R> = E Co 



+ J2 d aa u m ^ a ]nu + 0(u) m - a - 1 ]nw), 



s=0 



s=m+l 
s—m odd 



c as is undetermined if s > m and s — m is odd; 



(57) 



(58) 



\Q+m-s+l 



TVs m \ 

1 \2 2 > 



(s - m - 1)! (m - s + a)\ T(f - f - a) 



(if s > m and s — m is odd); 



(59) 



a-l 

E i 

i=o 

J even 



r(f + i: 



j!(a-l-j)! rU + l-a 



&OLS ~t~ 2 d c 



2 1 2 

if s < m or s — m is even. 



(60) 



Remark: Writing lntu in fl5T| ) instead of In ku, k some numerical constant, is arbitrary. 
In fact, in an application, k is likely to have physical dimensions (such as length). Changing 
k redefines the undetermined coefficient c as by adding a multiple of d as . 

2. Now we contrast the foregoing calculation with the parallel calculation of -R"/i from 
R^/i, assuming the latter to have the form fl5"7f). We use fl43f) and ( [44]) with k = 2: 



R$fi = 2 a R> + t J2, a {\ <t)R> da, 
Jo 



This time the j term vanishes for all j < |a — 1 (equivalently, j < [a/2]). Hence the 
contribution of small a will be of the form (as A — > oo) 

a-l 

i2= E ^A^'- 1 = 0(A- Q/2 ). (61) 
Hf] 

The integral of the remainder in ([h]) (over large a) is less than 

E K j x ~ j ~l a^-^ +j -hnada = terms of form (61) + 0(A^~^~5 in A), (62) 

r n J CTn 



3= f 



except that the second term is 0(X ™ 2 i (In A) 2 ) if the value j = |(a — m — 1) occurs. Here 
the relevant integral formulas are 

yV" 1 In Ac/A = A^GtTMnA-p- 2 ) + C ifp^O, / A^ 1 In Ac/A = i(lnA) 2 + C. (63) 
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Since m > 1, it follows that R and the similar terms in (|62|) are of at least as high order as 
the last term in (|62|). Hence there are, in effect, no undetermined constants of integration 
when the calculation is done in this direction (u — > A) . The remaining, explicit integrals are 



s=0 J(7 r 1 



1 



r(2j + 2) , 



1 



A a- 1 



+ E oda./ E ' ; A ' 1 

«=m4-1 z JfJ o . r„l 



j!(a-l-j)!r(2j + 2-a 

1 r(2j + 2) 



s=m+l 
s— m odd 



j!(a-j)!r(2j + 2-a) 



a 2 2 +J l n cr rfcr. 



Using fl6"3|), one sees that the values of s for which d as = yield terms proportional to 
c as A (m " s)/2 , but that when d as ^ there are, a priori, two types of term, proportional to 



c as A 2 2 + -d as X 2 2 ln A and a as A 2 2 . 



respectively. (Since 



m s . a ( a 

2 -2 +J> "2 + U" 1 



no logarithms are "created" by the integration as in (0).) Since it must be possible to 
recover the expansion (B^) in this way, there must be a numerical coincidence which causes all 
terms proportional to d as ln A to cancel. By the observation just made, this same numerical 
identity will cause all terms involving c as to cancel, if s — m is odd and positive. Thus these 
numbers c as do not affect at all the asymptotic behavior of R"fj, — as was to be expected 



from the fact that the latter does not determine them (pq) . Obviously, similar cancellations 
of logarithmic and c as terms must occur when is calculated from R^fj,, . . . , R^fi by 
Hormander's formula (|49|). 

The conclusion therefore is that if R®fi has the asymptotic behavior ([5j]), then R^fi has 
the behavior (|50|) with error term 0(A^~t~2 (In A) 2 ) in the worst case. The formulas for the 
coefficients are 



1 a— l 

- E 

2 4", 



-i 



-l + j 



T(2j + 2) 



j\(a-l-j)\ T(2j + 2-a) 

if s > m and s — m is odd; 



d, 



(64) 



a-l 



E(-i: 



a-j 



j\(a-l-j)\ 



T(2j + 2- a 
if s < m or s — m is even. 



2 a r 

as i * u as 



(65) 
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The conditions of consistency between the two calculations are (for < s < a, m > 1, 
a > 1) 



(from 



1 a-l im _ s _|_ ■ , i 



-2 



r(2j + 2) 



-i) 



(from (Eg) and (B3)); 



a-l 



3=LfJ 

a+m—s+l 



j\ (a - 1 
1 



r(2j 



a) 



r (! 



2 > 



1)! (m 



f -a) 



if s > m and s — m is odd 



r(2j + 2) 



x 



a-l 

E( 

i=o 

.j even 



\<x-3 



j\(a-l-j)\ T(2j + 2-a) 
s+j + 1)- 1 



+ 2 C 



■ 2 1 2 

if s < m or s — m is even 



and (H)); 



/ m 
a-j \ 2 



f+j + i)- 1 r(2 J + 2) 



+ 2° = 



3= t 



j!(a-l-j)! r(2j + 2-a) 

if s > m and s — mis odd 



(66) 



(67) 



(68) 



(the "numerical coincidence"). These can easily be verified for small values of a. 

In the Appendix it is shown that (|67D is valid for all values of s — m as a complex variable, 
except those where one of the factors is undefined. Then fl68D is obtained (for any a) in the 
limit as s — m approaches a pole of the second factor, and (^3) is obtained similarly from 
the derivative of fl67T). 



V. THE MEANS OF A STIELTJES INTEGRAL IN TERMS OF THE MEANS OF 

THE MEASURE 

The calculation in this section is the central lemma relating the asymptotics of Green 
functions at small t to the asymptotics of spectral measures of various kinds at large A. 
Let /i and / be functions of locally bounded variation, vanishing at 0, related as in (p9|) : 

/(A) = f a(a) dfi(a). 
Jo 

Assume that a(A) is a C°° function, and that it is well-behaved at the origin as described in 
the remark at the end of this section. Integrating by parts, one has 
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/(A) = a(AWA) - [ X a'(a)fi(cx) da. (69) 
Jo 

Continuing similarly, one can express d^ a f in terms of integrals whose integrands involve 
only <9^ Q /i and derivatives of a: 

ar/(A) = a(A)a A ->(A) + E(-i)^ T J 7— ^ jf (a - ay-^Xa^^) 



a+l rX raj-! .( 

£ / d<7 X ... / d<r(-lW 



a+1 



dia(a)d; a fi(a). (70) 



(The integral in the last version is j-fold.) This relation can be written 

01+1 ( n + 1 \ 1 r x 
R<*J(\) = a(\)R^(\) + \-<*J2(-iy( . J^-^j o (\-ay- l a a dia(a)Ry(a)da. 

(71) 

Proof: For a = 0, (|70|) is (|69|). Assume (|70|) for a and prove it for a + 1: First, 



dr~V(A)= [ X d^f(a)da 
Jo 



cr,_i 



A "+ 1 _|_ l\ r A 

a(a)a->(a)cia+ V(-1) J ( / da / da 1 ■ ■ ■ T ' dr dia(r)d; a /j(t) 

10 ~[ \ j J Jo Jo Jo 

The last term can be rewritten by (|33|) as 
Q+1 , /a + 1 \ 1 



a+1 / 

E(-i)T 

i=i V 



(A - aydia(a)d- a n(a) da. 



Now integrate by parts in both terms: 

0r-7(A) = o(A)^V(A) - fM^^-VWdff 

The last term equals 

E(-!)^ J J(— T)lJo (A - o-) j d a a(a)d~ a ~ 1 fi(a) da 

a + 2 fry _|_ l\ 1 i-A 

Using Pascal's triangle relation 

(T)+C-HT)< 
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one simplifies the expression to 

ar-V(A) = a(A)^-V(A) + E(-l) j (" t 2 ) //A - ^a(^W) ^ 

as was to be proved. 

Remark: In the foregoing it was tacitly assumed that <9^a(cr) was bounded as a { 0, so 
that there were no lower-endpoint contributions in the integrations by parts. Actually, since 

cT>(a) = 0{a a ) as a | (73) 

(as follows from (Q) and the boundedness of //), it suffices to assume that 

dy(a) = o(a' a ) as a | 0. (74) 



VI. HEAT KERNELS (LAPLACE TRANSFORMS) AND RIESZ MEANS 

The standard heat kernel of a second-order operator has a natural association with the 
variable we have called A, while the cylinder kernel (introduced in Sec. [TT]) is associated 
in the same way with the variable uj. The terms in the asymptotic expansions of these 
Green functions are in direct correspondence with those in the asymptotic expansions of the 
associated Riesz means. However, it is instructive to try to calculate each Green function in 
terms of the "wrong" variable, to observe how information gets lost, or needs to be resupplied, 
in passing from one quantity to another. Therefore, this section divides naturally into four 
parts. 

1. Given /i(A), a function of locally bounded variation on [0, oo), we consider 

roc 

K(t) = / e~ xt dn(\) (t > 0). (75) 
Jo 

As t I 0, we anticipate an expansion of the form 

oo 

K(t)~E^ rf+f ; (76) 



we shall demonstrate detailed equivalence of ([76]) with (|5D). 

In spectral theory, ( |75| ) has a number of possible interpretations. Let if be a positive, 
elliptic, second-order differential operator on an m-dimensional manifold. Then: (1) If the 
manifold is compact, K may be the integrated "trace" of the heat kernel of H, fi(X) being 
the number of eigenvalues less than or equal to A. (2) K may be the diagonal value of the 
heat kernel at a point x, fi(\) = E\(x,x) being the diagonal value of the spectral function 
(integral kernel of the spectral projection). In this case one knows that b s = if s is odd. 
(3) K may be the diagonal value of some spatial derivative of the heat kernel, [i being the 
corresponding derivative of the spectral function. (In this case, m in (^) depends on the 
order of the derivative as well as on the dimension.) If m — 1, fi is a Titchmarsh-Kodaira 
spectral measure. (4) K = K(t,x,y) may be the full heat kernel (off-diagonal), fi being the 
full spectral function. 
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Let us assume 



, ^ ^ — ^ m s _ , m a 



m Q 1 s 
2 



S=0 



and calculate K(t) from ([71]), with o(A) = e A ', 

d{a(\) = (-tye~ x \ 

and A — > oo in (|7T|). (See the first theorem of consistency, (|0).) In the present case all 
integrals in (|7ID converge as A — > oo, and hence the only term that survives in the limit is 
the one with j = a + 1, (A - a)^ 1 = X a + o(A°). So 



Kit) = (-1) Q+1 — / a a (~t) a+1 e~ (7t R a n(a) da 
a\ Jo 



a\ Jo 



m iOl 1 v 

2' 2 2 



s=0 



e~ at da. 



Using 



we have 



/■oo 

/ X'^e^dX = T(p)r p (p > 0), 
Jo 



K(t) = ± r{ ™J + " +1) a as rf+i + 0(r*+* 

s=0 r(a + i) 



m i a i 1 . 



(77) 



(78) 



(The undeterminable contribution from small cr is 0(t a+1 ), which is of higher order than the 
remainder.) Thus 



r(f -f + q + i) 
rr« + i) 



(79) 



holds for a > s, and formally for smaller a. (In the latter context a pole of r(|(m— s)+a+l) 
may be encountered. The significance of (|79|) then is that a as = 0. Thus, when m — s is 
even, a must be sufficiently large before a as will uniquely determine b s , even formally.) For 
a < s, (|79|) may be taken definition of au, . 

To check that the right-hand side of ( ff9|) is independent of a, write ( |50|) as 



1 a 
d\ a vW 



(80) 



and differentiate formally: 



1 



9r + VA)~^E(f-|+«)^Af-f 



By comparison with (RH) for a — 1, one has 
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1 / m s \ 1 



a! V 2 2 +a ) aas ~ (a-l)\ aa -^ si 



which may be variously rewritten as 



r(f-f + « + i) _r(f-§ + a ; 



r(« + i) I» 

(showing consistency of (|7|)), or as 



a a -i,s (81) 



1 / 771 S \ 

a a _ M = -( — -- + a) a as (a > 1). (82) 
a \ 2 2 J 

In particular, we can stem the profusion of constants by choosing a ss as the fiducial member 
of the family {a as }. We have 



r(f + | + i 



r( s + i) 



a,. 



and 



i> + i)r(f + f + i) 
r(f -f + a + i)r( s + i 



(83) 



a ss (84) 



(meaning 0, of course, when the first factor in the denominator has a pole). 

Remark: When m — s is even, a as is zero for small a, and the nonzero value for large a 
(hence the nonzero value of b s ) may be regarded as a constant of integration encountered in 
the passage from a low-order Riesz mean to a higher-order one. A case of particular interest 
is the local diagonal value of the heat kernel, for which b s = when s is odd but generally 
b s ytz o for s even. Here there is an essential difference between m even and m odd. If m 
is odd, then ao s is for odd s and nonzero for even s, and the constants of integration are 
also 0. If m is even, then a Qs = for all s > m, and the nonzero values of h s for s > m, 
s even, come entirely from constants of integration. This dimensional property of ao s is 



nitcg 

.EM 



reflected in the poles of the zeta function. 

2. It is of interest to consider how ( |76|) could be calculated from the cu-means, (0). We 
have 



<9j 




= e-» H { 


-2ut), 






= e~" H { 


\-2t + Auj 2 t 2 ), 


ft 




= e~ uH { 


[\2ujt 2 -8coH 3 ), 


ft 




= e~» H { 


[12t 2 -48w 2 t 3 + 16uH 4 ), 



(85) 



and in general the form 



If 

,2 



dl(e- u t ) = e-" H z^J- 2i i?-\ (86) 



i=0 



As before, when uj — > oo the formula (|71~D reduces to one term, 
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K{t)=YmxRZ / e-^d^r) 
where (|57j ) 



1 r°° o 
al Jo 



+ £ d as oo m - s \nuj + 0(uj m ~ a ~ l \nu). 



s=0 



s=m+l 
s—m odd 



The relevant integrals this time are (p > 0) 



(87) 



e -u> t w 2 P -i lnujduJ = Y(p)hp(j>) - \nt]t- 
o 4 



where 



satisfies 



if>(p) = d p lnT(p) 



i>(p + 1) = i>(p) + - 
p 



1,2,...), 



ip[- + n 



"7 



2k 



(n = 0,l, 



(88) 
(89) 
(90) 
(91) 
(92) 

7 = +0.57721566490 (Euler's constant). (93) 
From (RBI) we see that the contribution of small oj is 

Q(j.[(a+2)/2}^ and the i ntegra i Q f tne 

error term out to to = oo is 0(t~ T+f+2 hit); the latter is dominant, as in (ff8|). Therefore, 
the remaining, explicit terms should reproduce exactly the summation in ([7^). This requires 
cancellation of all In t terms; comparison of (^) with (|88|) shows that this will entail cancel- 
lation of all c as terms with s — m odd and positive, just as in the calculation at the end of 
Sec. [IV] (and also, coincidentally, of all terms containing 7 + 2 In 2). We illustrate by working 
out the case m = 1, a — 3, using the last of Eqs. fl35l): 

Kit) ~ - / e~ w \l2t 2 - 48cu 2 t 3 + 16oA 4 ) 
6 Jo 

x [c 30 lu 4 + c 3 ilu 3 + c 32 uj 2 + c 33 uj + d 32 uo 2 \nu]duj 

= c^r 1 / 2 [r (I) - 4r (z) + f r (§)] + c 31 [r(2) - 4T(3) + |r(4); 
+ C32 t 1/2 [r (I) - 4r (I) + fr (f )] + c 33 t [r(i) - 4r(2) + |r(3) 

- id 32 (7 + 2 in 2 + \nt)t^ 2 [r (I) - 4r (|) + §r (I) 

+ 4 2 t 1/2 



n— 1 j 
fc=l K 



21n2 + 2j] 

fc=i 



r(| 



3 1 V2 y '^45 ± V 2 



2r (I) C30 r 1/2 + c 31 + f r (§) d 32 * 1/2 - |c 33 t. 
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Moreover, the numerical coefficients agree with those computed from (|78|)-(|79D and fl64|)- 
3. Next consider the quantity 

POO 

T(t) = / e'^dfj, (t>0), (94) 



where \i = //(A) = fi(u; 2 ) = K{t) bears the same relation to the heat kernel of the 

operator H which T(t) bears to the kernel of the operator exp(— \f~H t), the cylinder kernel. 
That operator solves the elliptic partial differential equation 

d?ip(t, x) - Hip(t, x) = Q (t > 0) (95) 

in a cylindrical manifold of dimension m + 1, with inhomogeneous Dirichlet data on the 
m-dimensional boundary surface t — and a decay condition as t — > oo. 

The calculation of the small-t expansion of T(t) from the cu-mean expansion (57) starts 
off in precise analogy to the previous calculation of the expansion of K(t) from the A-mean 
expansion flSOf). One obtains from ([71]), with u; in the role of A, 

1 POO 

T(t) = (-l) a+1 — / T a (-t) a+1 e- Tt R«jldT 
a\ Jo 

f2c as u m - s+a + J2 d as u m - s+a \nu + 0{u m - l \nu) 



t a+l}_ 

a\ Jo 



s=0 s=m+l 
s— m odd 



-LOt 



du. 



In addition to (fTfi) in the form 



POO 

/ t^-^-^du; = T(p)r p 
Jo 



we need a transformation of (1881) 



/ e- ut u p - 1 hiudu= r(p)[V»(p) -lnt]t~ p (p > 0). (96) 

*/ 

The integral of the error term is of order 0(t _m+a+1 hit), and the contribution from small 
u is 0(t a+1 ), hence of higher order than the error term. Evaluating the integrals of the 
summations, we get 

T(t) = ± n™- s + <* + !) +s 

U t r( a + i) 

+ E ~ * ± " + 1} mm - 8 + a + 1) - hit] rf a .t-"+' 

s=m+i l (a + 1J 

s— m odd 

+ 0(r m+a+1 lnt). (97) 
Thus, if we define notation by 

oo oo 

T(t)^Y / e s t- m+s + £ / s t- m+s lnt, (98) 

s=0 s=m+l 
s— m odd 
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we have 

T(m — s + a + 1 



if s — m is even or negative, and 

r(m — s + a + 1) 



r(a + i} 



T(a + 1) ~ us (99) 



[c as + ip(m - s + a + l)d as ], (100) 



/. = - r(m ~ 8 + ° + 1) d as (101) 
J r(a + i) v ; 

for s — m odd and positive. These equations are rigorously valid for a > s, and hold formally 
for smaller a — that is, they can be used to define c as and d as for a < s, with one exception: 
When 

m-s + a<0, (102) 

the T function in the numerator has a pole. In the context of ( |99|) or ( |101| ), this is understood 
to force c as = or d as = 0. However, one may not conclude that the c as in ( |100] ) is zero in 
this situation, since ip also has a pole; indeed, we know that c as is generally nonzero in the 
case a = 0, m odd, where the relation of R^fi to R^fi and the heat kernel K(t) is trivial. 
The right-hand sides of (p9|) (p70T|) must be independent of a. To verify this, rewrite 



as 



dZ^~-,ilc as u m - s+a + - J2 d as Lu m - s+a \nu (103) 



a - s=0 a - s=m+l 

s— m odd 

and differentiate: 



+ — dascu™- 8 ^- 1 + — (m-s + a)d M w m -» +tt - 1 lnu;. 



a - s>m+l a - s>m+l 

s— m odd m odd 

This yields the recursion relations 



1 1 

— -Am - s + a)c as = rrTC Q ^i iS (104) 

a\ [a — 1)1 



for s — m even or negative; 



1 11 

— -(m - s + a)c as + — -d as = — - c a _i jS , (105) 

a\ al [a — 1)1 



— (m-s + a)d as = - — l —— d a _ 1)S (106) 
al [a — 1)1 
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for s — m odd and positive. Multiplying by T(m — s + a), we immediately notice consistency 
of (ELOD and fllO§ with © and (|MD [cf. ©], while (|105|) becomes 

T(m — s + a + 1) T(m — s + a) Y{m — s + a) 

: 77 c as H 777 j 77 "as — 777 n c a~l,s • 

r(a + l) T(a + 1) r(a) 

This, with ( |90| ) and ( |106| ), implies 



r(m-s + a + l) r ,. , T(m-s + a), . , . 

fTT — -~r\ [c as + ip{m - s + a + l)d as \ = — — [c a _i, s + ip(m - s + a)d a - 1:S \ 

1 [a + I) 1 [a) 

establishing consistency of (|100| ). 



The form of the recursion relations analogous to (|82|) is (for a > 1 



c a -i )S = -(m- s + a)c as (107) 
a 



if s — m is even or negative; 



c a -i, s = —{m - s + a)c as + -d as , (108) 
a a 



d a -i, s = -(m- s + a)d as (109) 
a 

if s — m is odd and positive. Note the consistency with the remarks made in connection with 
( |102| ): If s — m is odd and positive, then c as for a value of a satisfying ( |102| ) will generally 



be nonzero and proportional to d as for values of a violating ( |102| ). Otherwise, c as and d D 
will vanish when ( |102| ) is satisfied, because the factor (m — s + a) in ( |107| ) or ( |109| ) will have 
vanished for some larger value of a. 

Let us express all the coefficients in terms of c ss and d ss : In analogy to (|8"3"|), we have 



r( s + i) 

if s — m is even or negative; 



Css (HO) 



r( s + i, 

if s — m is odd and positive. The analogues of ([84]) are 



r ^ +1) [c, + ^ + iK s ], (in) 
rs + i 



r(m+l) , , s 

/. = , ' dss (H2) 



r(m-s + a + l)r(s + l) 
if s — m is even or negative; and in the contrary case, 



T(a + l)T(m + l) 

c ss (H3) 
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r(a + l)r(m + 1) 
T(m-s + a + l)T(s + l 



-{c ss + [ip(m + 1) - i){m - s + a + i))d ss }, 



(114) 



(In 



r(a + l)r(m + 1) 
T(m-s + a + l)T(s + l) 



(115) 



If m — s + a > 0, then ( [114j) is unambiguous, and the difference of ip functions is calculable 
from ( p0| ) or (pT|). But when ( |102|) holds, extra work is needed to obtain a usable formula. 
In that case iteration of (|108|) leads to 



■I) 1 



- a -i£(£ - m - a ) r (« + i)r(m + 1) 
r(s + i) 



(116) 



for a < s — m, s — m odd. This formula also can be obtained from ( |114[ ) by comparing 
residues of T and ^ 



^>(e — n) r'(e — n) 
rfe-n) = r(e-n) 2 



\n-l 



n! + 0(e) 



(n = 0,l, 



(117) 



(see (|§). 

4. Finally, we investigate the relation between T(t) and the A-means. For use in (|71~D 
we find 



d x (e 



d*(e- txl/2 



e -ov a (-fa- 1 / 2 

e" tAl/2 (i^A- 1 + itA- 3 / 2 ) , 
a 3 ( e -' Al/2 ) = e - tAl/2 (-§i 3 A 3 / 2 - |t 2 A" 2 - |tA~ 5 / 2 ) , 



<9?(e 



(118) 



and in general the form 



9((e 



-^A 1 /2^ 



-a 1 / 2 



1=1 



(119) 



Note that ( |74|) is satisfied. As in the other cases, there is only one relevant term in (|Tl): 



T(t) = (-l)^ 1 ^ [°° d« +1 (e- txl/2 )\ a miid\ 
a'. Jo 



-1 



a! jo 



js=0 



i+a 



+ 0(A 



m j_q 1_ , 

2^2 2 



dA. 



(120) 



However, unlike the other three cases, this time the unknown contribution from small A is 
0(t) — its order does not increase with a. Furthermore, the integrals of the individual terms 
in ( |120|) do not converge at the lower limit, although the integral as a whole is convergent. 
Consequently, we must deal with integrals of the form 



X~ p d\ = 2 / e- tw u- 2p+1 duj. 

U) 



(121) 
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If p < 1, we may set lu = and use (|77|) ; if p > 1, we may use 

e-^" n cL; = ul~ n / e-" oiu M- n du, (122) 



/oo 
e- tw u- n duo = E n (t), (123) 



= _ 7 _ i n t _ £ L_ Lip, (124a) 



^(*) = M T whlnt + V'(n)]- £ (124b) 



(n-1)! 1 f-T, 0'-n + l)j! 



It can be seen that the integrals give rise to terms of the form fl98|), but without the restriction 
that s—m be odd in the logarithmic terms. Moreover, the contribution from A < 1 is analytic, 
hence can be expanded as a series of positive integral powers of t with unknown coefficients. 
From ( |58[) and (j99|)- (|100|) , we would expect the coefficient of t~ m+s to be determinable 
if s — m is even but not if s — m is odd. Therefore, we anticipate conspiracies among 
the coefficients in ( |119| ) which will eliminate both the logarithms and the undeterminable 
coefficients when s — m is even. Let us verify this for the case m = 1, a = 3: 

Tit) = - r e- txl/2 ( -t A \~ 2 + -t 3 A" 5 / 2 + -t 2 \- 3 + -t\~ 7 / 2 ) 
w 6 Jo V16 8 16 16 J 

x [a 30 A 7/2 + a 31 A 3 + a 32 A 5 / 2 + a 33 A 2 + 0(A 3 / 2 )] d\ 

1 f°° 

= -JE / e- tuJ (a 30 t A u A + 6a 30 t 3 u 3 + a 3l t A uo 3 + 15a 30 tV + 6a 31 t 3 cu 2 + a 32 t A uo 2 
48 jo 

+ 15a 30 tuj + 15a 3 it 2 u; + 6a 32 t 3 a; + a 33 t 4 u + 15a 31 t + 15a 32 t 2 + 6a 33 t 3 
+ I5a 32 tu~ l + lha 33 t 2 u~ l + 15a 33 tcr 2 + 0(iuT 3 )) du 



-^{l05a 3 o* _1 + 48a 31 + 23a 32 t + 7a 33 t 2 



+ 1 e~ tw {l^t 2 uo- 6 )[0{uo & )-a 32 uo b }ujduo+ e-^ilhtu^Oiu^oj du + 0(t 
Jo Jo 

+ 15o32tSi(woO + lba^E^uot) + l5a 33 u^HE 2 (uj t) + / e- tuJ 0(tu~ s )duj 

Jwn 



(Here the first four terms come from integrating all terms in the middle member involving uj 
to a nonnegative power. To represent accurately the integral over small u of the remaining 
terms, it has been necessary to go back to the first member and to note that the quantity 
in square brackets there is 0{uj & ) by virtue of ( |1 20|) and (ff3|). Of the terms involving 



t\~ 7 / 2 , the first two have already been accounted for in the explicit integrations, and the 
remainder of the square bracket is still 0(u 6 ) as u j 0; of the terms involving t 2 A~ 3 , the 
first three have been accounted for, and the subtraction of the third of these needs to be 
represented explicitly in our formula; the terms associated with the rest of d A (e~ txl/2 ) are 
0{t 3 ).) Substituting from (|122|) -( |l"2"4T ), and noting that the two quantities represented by 
"0(u; 6 )" in the foregoing expression are the same, we obtain 
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T(t) = ^|l05a 3 ot~ 1 + 48a 3 i + 23a 32 t + 7a 33 t 2 - I5t 2 a 32 u 

+ / 15tO(u )[ujt + 0{t 2 ) + l-ut]dcj + 0{f) + tO{u Q 2 E 3 {u} t)) 

J 

+ 15a 32 t[-7 - In wot + co t + 0(t 2 )) + 15a 33 t 2 [-7 - hiu t + 0(t)] 
+ 15a 33 cJ 1 t[-cj t(-lncj t + ^(2)) + 1 + 0(t 2 )]J 

= -^|l05a 30 t _1 +48a 3i - 15a 32 tlnu; t 

+ t 23a 32 + 15 / 0(uj°)duj - lh^a Z2 + Iha^ 1 
Jo 

+ t 2 [7a 33 - 15a 32 ^ + 15a 32 w - 157^33 + 157^33 - 15a 3 3] 
+ 0(f) + u 2 t [O (I - u t) + 0(t 2 hit)] | 

= — {lOSa^™ 1 + 48a 31 — 15a 32 tlnct; )f: + const, x t — 8a 33 t 2 
48 

+ 0{f hit) + 0{uo 2 t) + 0{uqH 2 )}. 

The bothersome term at the end could be removed by taking ujq —>■ oo. Now this term and 
its predecessor can be traced to the first two terms in a Taylor expansion of e~ tuJ in (|1 23|) 



with n = 3. In analogy to the argument leading from Q55D to (p7D, we may transfer this part 
of the integral from the high-cj account to the low-cu account; this simply amounts to adding 
something to the "0(o; 6 )" term, and the same cancellation of 0(t 2 ) terms just observed for 
that term will occur for this other contribution as well. Thus we have derived the expected 
form, 

T(t) = ffaaor 1 + a 3 i - ^a 32 tlnt + e 2 t - |a 33 t 2 + 0{t 3 \nt), 

where e 2 is undeterminable from the asymptotic expansion of R\[i. Furthermore, the coef- 
ficients agree with those calculated from (pgp- flnnp and (B3)-(Bn) 



Summary: The coefficients b s in the asymptotic expansion of the heat kernel, K(t), 
are in one-to-one correspondence with the "diagonal" coefficients a ss in the asymptotic 
expansions of the Riesz means with respect to A, -R^yU. The nondiagonal coefficients, a as , 
in the expansions of the R"/j, differ only by numerical factors from the a ss ; because these 
factors may vanish, it is not possible to express a ss in terms of a as if a is too small. Similarly, 
the coefficients e s and f s in the expansion of the cylinder kernel, T(t), are in one-to-one 
correspondence with the diagonal coefficients c ss and d ss in the expansions of the Riesz 
means with respect to u, -Rfj/i. The connection in this case involves a two-termed equation, 
( p. Again the c as and d as can be expressed in terms of the c ss and d ss , but not conversely 



if a is too small. Finally, and perhaps most significantly, the c ss and d ss (or the e s and f s ) 
contain information which is not contained in the a ss (or the b s ) (but not conversely). The 
c ss for s — m odd and positive are "new spectral invariants" independent of the a ss . In the 
case of the heat and cylinder kernels, the c ss contain nonlocal geometrical information, while 
the a ss are strictly local. 
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VII. CONCLUSION 



Motivated by applications in spectral asymptotics and quantum field theory, we have 
investigated Riesz means in this setting: Two functions, K(t) and T(t), are related to a 
function //(A) by the Stieltjes integrals 

/•oo roo 

Kit) = / e~ xt dn, T(t) = / e-"* d/i, (125) 
Jo Jo 

where u = vA, and they have the asymptotic expansions 

oo oo oo 

K(0~EM~ f+f , TW~^e s r +s + E fst~ m+s lnt, (126) 

s=0 s=0 s=m+l 

s—m odd 

as t 1 0, for some m G 7L + . The behavior of as A — ► +oo is characterized by the constants 
a ss , c ss , 4s , where 

#>(A)~X>- Af ~ f > i?>~Ec QsW m - s + E ^ m ' s ln^, (127) 

s=0 s=0 s=m+l 

s—m odd 

R"fi and -R^/ 1 being the Riesz means with respect to A and u, defined by, for example, 

Ry(u 2 )= f (l-lV^(r 2 ) (128a) 

Jt=0 V LU/ 

= a\ Lu~ a r dn ■ ■ ■ H' 1 dr a /i(r Q 2 ). (128b) 



o 

We have shown that a ss and b s contain the same information, being related by fl83|) , and 
that (c ss ,d ss ) and (e s , f s ) contain the same information, being related by ( p. 10| )— ( pTT2|) . The 
inverses of the cited formulas are 

F(s + 1) K; (129) 



r(f + | + i 



d ss = — -=-r — — — "T - f s if s — m is odd and positive; (130) 
T(m +1) 



c ss = — — -r[e s + ip(m + l)/ s ] if s — m is odd and positive, (131) 

1 \ Tfl ~\- 1 J 



e s it s — m is even or negative. (132) 



r(m + l) 



Finally, these two collections of quantities are related to each other, but in an asymmet- 
rical way. Specializing d58|)-(|60|) and (|64|) — (|65D , we have 
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c ss is undetermined by a ss if s — m is odd and positive; 



(133) 



d. 



-1) 



m+l 



i r(f-f) 

(s-m-l)!m! T( - f - f 



■ a ss if s — m is odd and positive; (134) 



+ 



s-1 

E 

j=0 
j even 



("I) 



^ (m-a + j + l)- 1 T(f + |) 

3 \(s-i-j)\ r(| + |- s ) 



if s — m is even or negative; 



(135) 



1 

•- E 

2 4", 



. (™ _ £ + ? ' + ;p 



r(2j 



T(2j + 2- S ) 
if s — m is odd and positive; 



(136) 



2 S + E 



. fs _ £ + j + i 



,-1 



r(2j + 2) 



i> 1 ./)! T(2j + 2- s ) 

if s — m is even or negative. 



(137) 



These relations obviously induce similar relations between b s and (e s , f s ). The quantities c ss 
(s — m odd and positive) summarize the information contained in the asymptotics of T but 
not in that of K. 

In the context of spectral theory of differential operators, the quantities b s , e s , f s are 
more accessible to calculation, because K(t) and T(t) are Green functions or their traces, 
but a ss , c ss , d ss are more fundamental, since fi(\) represents a spectral decomposition (the 
integral kernel of the spectral projections in the local case, the density of eigenvalues in the 
traced case). We expect that further work will bring the technically more difficult case of the 
Wightman function (and other Green functions of wave equations) into the same picture. 

This picture clarifies and unifies results that have appeared in the literature piecemeal 
or imprecisely. One sees that the small-t expansion of K is not related in a one-to-one 
way with the formal large-A expansion of /i (i.e., the ao s ); rather, the heat coefficients b s 
can involve constants of integration from the Riesz averaging of /i over small A. This ex- 
plains how the formal expansion of fi (and the associated zeta-function poles) can be strik- 
ingly different in even and odd dimensions, while the heat-kernel expansion is notoriously 
dimension-independent. (See Remark below (|84]).) Further constants of integration appear 
in the passage from A-means to w-means; these carry nonlocal geometrical information in 
the spectral context. They are intimately related to the potential appearance of logarithmic 
terms (f s , d ss ). On the other hand, the numerical relationships among the various series are 
independent of the spectral application; in particular, they will be universal for all (positive, 
elliptic) operators of a given order and dimension. (Only order 2 has been treated here.) 

To conclude, we apply some of these formulas to the simple examples studied in Sec. |TJ 
and compare the results to the known spectral densities for those cases. In these examples, 
m — 1. 
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Case M = R. 
Formula flTU| ) for the heat kernel states that 

b s = (47r)- 1 / 2 5 s o if V = x, (138a) 

b s = for alls ify^x. (138b) 
From (129), then, one has for the Riesz means of the spectral function 

- 5 s0 ify = x, (139a) 



"ss 



7T 



a ss = for all s if y ^ x. (139b) 

The eigenfunction expansion in this case is the Fourier transform, so the exact spectral 
density is 



dE x (x,y) = — Y'e ik ^d\k\ = -cos[uj(x-y)}duj (V\ = uj=\k\). (140) 

sgn k 

Thus 



E x (x,x) = - = y^, (141a) 

7T 7T 



_, , . 1 sm[u(x - y)} . . 

E\(x, y) = ^ ^ lfy^x. 141b 

7r x — y 



(Constants of integration are fixed by (|28|) .) Clearly (|141a ) is consistent with ( |139a|) . Equa- 



tion ( |139b| ) indicates that ( |141 b| ) is "distributionally small", or decays rapidly at infinity in 
the Cesaro sense.00'i To see this by elementary, classical methods, one would use (j36[), 

l — s 

and verify that this object is o(X~) as A — > oo. This can be done by repeated integration 
by parts (most easily after changing variable from a to r = \fo). 
For the less familiar cylinder kernel, ([121) gives 



/. = 0, (143) 



e s = - 5 s0 ify = x, (144a) 

7T 



and if y ^ x, 
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7T (x - y)' 



if s = or s is odd, 

if s is even and positive. 



From (|130|) - (|132|) , this corresponds to the cj-means 



0. 



(144b) 



(145) 



1 A 

- o s0 

7T 



if y 



x. 



and if y ^ x, 







(_l)l+i 



if s = or s is odd, 
if s is even and positive. 



(146a) 



(146b) 

7r (x-yY 

One observes that (|139|) and (|145|) - (|146|) are related precisely as prescribed in (|133|) - (|137|) . 
Again, ( |145| )-( |i~46| ) could be established directly from ( |128a| ) and ( |141| ) by integration by 
parts. (The earnest student who actually attempts these tedious exercises will find that 
certain endpoint terms that vanished in the previous case will produce the nonvanishing c ss 
in this case.) The present case, however, is more easily treated by ( |128b| ). 



Case M = R H 



Let us consider only the Dirichlet boundary condition. The heat-kernel asymptotics, and 
hence the A-means a ss , are again trivial. The effect of the boundary is seen (as dictated by 
( |133|) ) only in the even-order coefficients of the cylinder kernel and the associated tu-means: 
from (|T5l), for s > 0, 

(-l)»/2 1 



vr(2x) 



c ss if y = x and s is even, 



(_l)«/2 



7T 



(a; - yY (x + y) 1 







— c ss if y 7^ x and s is even, 
si 



if s is odd. 



The eigenfunction expansion in this case is the Fourier sine transform, so 

2 

dE\(x, y) = — sin(a;a;) sin(ujy) du. 

71 

Let us check consistency only for the diagonal values. We have 



Ex(x,x) = — 



sin (2ujx) 
2ttx 



(147a) 

(147b) 
(148) 

(149) 
(150) 



The first term is the same as ( p. 4 lap and need not be discussed further. The A-means of the 
second term are shown to vanish exactly as for ( |142j ). The w-means of the second term can 
be calculated by (|128b| ) exactly as for ( |146b|) , and they reproduce (|147a|) . 
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Case M = S l 



The eigenfunction expansion is the full Fourier series, so 



dE x (x,y) 



2L 



sgn k 



Jk(x-y) 



n=0 



E* 1*1- 



nn 



d\k\ 



duj. 



;i5i) 



The traced expansions of the Green functions correspond to the eigenvalue distribution 
function, 



dN(X) 



dE\(x, x) dx 



-L 



du. 



(152) 



The calculation of Riesz means leads to sums that can be regarded as trapezoidal-rule 
approximations to integrals of the types already considered. Therefore, the new features of 
this case can be sought in the Euler-Maclaurin formula (Ref. Sec. 1.8) for the difference 
between the sum and the integral. The terms in that formula involve only the odd-order 
derivatives of the integrand at the endpoints of the interval. In an integral such as (|128a|) 
or (142), the derivatives of the factor such as (1 — t/uj) s to all relevant orders will vanish at 
the upper limit, so only the lower endpoint can contribute. In the A-means (|142|) , where this 
Riesz factor is (1 — r 2 /A) s , the odd-order derivatives all vanish at by virtue of factors of 
r or sin[r(x — y)\. This is necessary for consistency with the trivial heat-kernel expansion, 
(PH) and fl2Tf). For the cu-means, we need to calculate 

p d v \ ( rV 

cos [t(x — y)} 



1 

L 



drP 



l- T - 



(153) 



When y = x, this derivative is of order u> p ; in fact, it equals 

(-l)Vs! 
Lp +1 ujp{s -p)\ ' 



(154) 



(In the off-diagonal case, the Euler-Maclaurin series does not yield an expansion of the 
desired type — at least, not without a resummation, which we shall not attempt here.) 
To read off c ss (for s > 2) we need to look at p = s — 1 (see (|127| )) and, according to 
Euler-Maclaurin, multiply by (— l) p B s /s\, where B s is a Bernoulli number. Thus we arrive 
at 



TT 



3-1 
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s-1 



L s 1 a L s s\ 
This is, in fact, the same as (|TTD-(|i"8"D with y = x, because 

1 sinh,?: z 1 
= coth - — - 

2 cosh z — 1 2 z 



oo r> 



s=2 



SI 



(155) 



(156) 



Be = if s is odd and s > 2. 



(157) 
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Case M = I 



For Dirichlet boundary conditions the eigenfunction expansion is the Fourier sine series, 

so 

2 00 / nir\ 
dE\(x, y) = — sin (lox) sin (my) 5 ( uj — ) du (158) 

L n =i \ L 

and 



00 / rar\ 

dN(\) = Y,8(u- — ) du. (159) 

n=l ^ L ' 



The novel feature of this case is that the trace (|26|) of the heat kernel contains a nontrivial 
term, of order t°, and a similar term appears in the trace fl2~?|) of the cylinder kernel. The 
corresponding Riesz means are 

an = h = - \ , en = ex = - \ . (160) 

In an Euler-Maclaurin calculation these terms arise from the innocent fact that the initial 
term of the trapezoidal rule, 

|5(w-0)dw, (161) 

is missing from ( |159| ) and needs to be subtracted "by hand". 

The reader will have noted that in most of these instances it is easier to obtain the 
Riesz means from the Green functions via ( 129 )-( |132|) and Section |lj than to calculate the 



Riesz means directly, even though the spectral functions are known exactly. (Since the 
direct calculations are merely consistency checks for us, we have carried out some of the 
calculations only far enough to demonstrate the matching of the first few terms, rather than 
constructing complete proofs. The claims in this section are to be understood in that spirit.) 
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APPENDIX: PROOF OF (|67|) 

BY R. A. GUSTAFSON 



We shall prove the following identity: 



lr- 



X 



^ l ) j\(a-l-j)\r(2j + 2-a)^ 



,V J j!(a-l-i)!r(i(i + l)-a) + 

j even 



(Al) 



where a > 1 is an integer and z GC is not a pole of either factor. The key fact used is the 
3_F 2 transformation formulaEl, 



3 r 2 



( a, b, —n 



V e,/ 



(e - a) w (/ - a) T 
(e)n(f)n 



/l- 



3^2 



s, a, — n 



\ 1 4- a — / — n, 1 + a — e — n 



(A2) 



where s = e + / — a — b + n. 

First we compute the sum in the first factor of (|A1|). With a little work one finds 



a-l 



V ( 1)a - 3 - (^ + J + 1 )" 1 r(2j + 2) 



j! (a-l-j)!r(2j + 2-a) 



;-l)(2«-l)! 



3-^2 



Case 1: a is even. In 
expression (|A3|) becomes 



(f + a) [(a- l)!] 2 

set a 



.« _ a —°L _|_ I _« _|_ 1 
2 ' 2^2' 2 ^ 



(A3) 



a + 1, —a + 



2 

■~ + |, 6 = — | — a, and — n = — ~ + 1; then the 



(-1)(2«- 1)! 



2 2 ~ 2^a/2-lv 2 /a/2-1 



1 4- - -- + 1 

' 2^2' 2 ^ 



(f + a) [(a - l)!] 2 ( - f - a + -a + |) a/2 _ x 



3^2 



- 4- - 2 

2 ^ 2' Z 



1 . (A4) 



Now observe that 



3-^2 



/ 1 _« _i_ I -« + 1 

l ' 2 ~2' 2 ' 

* i 3 



VI + 1,2 



a _ 1 

2 2 



)( 



Z + 



a 1 \ / a\ 

2 2/V II 



3 r 2 



I i a 1 a 

' l i 2 2' 2 



V - + - 1 

V 2 ~ 2 ' 

a 1 a 

'2 2' 2 



;i -i 



2^2 



(A5) 



Hence the expression in flA"3]) becomes 
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(-l)(2a-l)! 



>+§) 



+ a)[(a-l)!]»(-f-|)(-f) 



1 

2> 



V - + - 

\ 2 ~ 2 



X 



2)a/2-l( 



2 /a/2-1 



( 



« + 1 )a/2-l(-« + 



2)0/2-1 



(A6) 



By the Gauss summation theorem for 2 i*i , this is 



(-l)(2a-l)!(z+|)(- 



2 T 2/a/2-l l 



2 ) a/2-1 



(f + a) [(a - 1)!] 2 ( - f - §)( - f )( - f - a+ 1)^ - a+ 1) 



a/2-1 



X 



r(§ + |)r(§ + a + i) 



r(f + f + i)r(| + f + i) 

After carrying out a long, easy calculation we obtain: 



a-l 

E(-i) c 

3=0 



(iz + j-iyi r(2 J + 2) 



V2 



v2 2) a/2 



j!(a-l-j)!r(2j + 2 



a: 



+ f + l)a 



/2 



recalling that a is assumed to be even. 

Case 2: a is odd. We use (|A3|) , now setting a — ,, ...... 

Doing a calculation very similar to the previous case, we obtain 



■f + 1,6 



^ 1V ,_, (^ + i + i)' 1 r(2i + 2) 



2 a 



+1) 



2^a/2+i 



2 1 o 1 1_\ 
2 " r 2 " r 2 /a/2+1 



a and — n 



- 2 Q . 



(A7) 



(A8) 



- + - 

2^2' 



l-i)!r(2j + 2-a 

We now turn to the sum in the second factor of ( |A1| ). With some work one obtains 

(z+j + i)- 1 r(i( J + i)) 



(A9) 



Q-l 



j even 



(z + l)(a-l) 



a 771 
3-^2 



2 + 1 

2 ~ 2 : 



2' 



.2 + 1 
2^2' 



- + 1 

2 ^ 



;i 



-a 



Case 1: a is even. Again we apply (|A2|) , setting a = — 
After a computation quite similar to the above we obtain 



a-l 

y (-i) a ' j 



( g + j + i)-ir(i(j + i)) 
_-j)\ra(j + i)-a) 



- + - b 

2 ~ 2' 



:i + f + 1) 



(A10) 

+ i and -n = -# + 1. 



a/2 



- + - 



- 2" 



2 /a/2 



7 even 



Case a zs odd. Set a = — | + 1, 6 = | + | and — n = — | + | in 



(All) 



After another 



computation very similar to the previous ones we find 

a-l 

E (-!) a ~ j 



(^+i + i)- 1 r(l(i + i)) 



2~a^2 + 2 + 2)a/2+i 



j! (a - 1 - j)\ r(|(j + 1) - a) ~ " (f + |) Q/2+| 
Substituting (|A8|) - (|A12|) back into (|A1|) yields the desired identity 



- T 



(A12) 



j even 
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